Can binary mergers produce maximally spinning black holes? 
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Gravitational waves carry away both energy and angular momentum as binary black holes inspiral 
and merge. The relative efficiency with which they are radiated determines whether the final black 
hole of mass Mf and spin Sf saturates the Kerr limit (x/ = Sf/M'j < 1). Extrapolating from the 
test-particle limit, we propose expressions for Sf and Mf for mergers with initial spins aligned or 
anti-aligned with the orbital angular momentum. We predict the the final spin at plunge for equal- 
mass non-spinning binaries to better than 1%, and that equal-mass maximally spinning aligned 
mergers lead to nearly maximally spinning final black holes (x/ — 0.9988). We also find black holes 
can always be spun up by aligned mergers provided the mass ratio is small enough. 



Recent breakthroughs in numerical relativity have led 
to successful simulations of the inspiraljinergcr, and ring- 
down of binary black holes (BBHs) [l|, l2, l3|. Despite 
this progress, approximations underlying these simula- 
tions limit them to mergers with large but submaxi- 
mal initial spins {xi ^ 0.9) and order-unity mass ra- 
tios {q = railm\ > 1/6). Yet the maximally spinning 
regime is of considerable theoretical and observational 
interest. Maximally spinning black holes barely manage 
to hide their singularities within their event horizons, 
and observable, "naked" singularities are expected for 
spins X > 1- While infinitesimal processes like steady 
accretion cannot produce naked singularities [j|, Pen- 
rose's cosmic-censorship conjecture [5| that such singu- 
larities can never be created has not been proven for 
comparable-mass BBH mergers. Black-hole spins are not 
a purely theoretical concern; they can be measured by re- 
verberation mapping of iron Ka fluorescence in the spec- 
tra of active galactic nuclei (AGN) [6j. This technique 
has been applied to XMM-Newton observations of the 
Seyfert 1.2 galaxy MCG-06-30-15, leading to a measured 
spin X = 0-989tooo2 very near the maximal limit Q. The 
spins of supermassive black holes (SBHs) also offer im- 
portant insights into their formation. Some expect that 
highly spinning SBHs will only be found in gas-rich sys- 
tems like spirals [8[. Others suggest that gas accretion 
occurs through a scries of chaotically oriented episodes, 
leading to moderate spins (x ^ 0.1—0.3) lower than those 
expected from comparable-mass mergers in gas-poor el- 
lipticals [9|. Comparing measured spins in spirals and 
ellipticals will distinguish between these two scenarios. 



In the absence of reliable simulations of maximally 
spinning BBH mergers, we must rely on various approxi- 
mations in this imp ortant regime. Hughes and Blandford 
(hereafter HB) [lO| assumed that the energy and angular 
momentum radiated during the inspiral stage dominates 
that radiated during the comparatively brief merger and 
ringdown. They then used conservation of energy and 
angular momentum to equate these quantities when the 
BBHs reach their innermost stable circular orbit (ISCO) 



to the mass and spin of the final black hole. 
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Here -B(x) is the energy per unit mass of a test particle 
on an equatorial orbit of a Kerr black hole with spin pa- 
rameter X ^i^d Lorb (x) is the corresponding dimensionless 
orbital angular momentum. While E{x) and Lorb(x) are 
only available in the test-particle limit (m-i <^ m\) [ll| . 
HB extended these exact, analytic results to comparable- 
mass mergers. While this approach reproduces the test- 
particle limit by design, it is not symmetric under the 
exchange of labels "1" and "2" and leads to supermaxi- 
mal spins (x/ > 1) for near maximal mergers. 

Buonanno, Lehner, and Kidder [l^l (hereafter BKL) 
remedied these problems by (1) conserving mass, (2) in- 
cluding the initial spin of the smaller black hole, and (3) 
using the spin parameter of the final black hole X! rather 
than that of the more massive of the initial BBHs xi to 
determine the orbital angular momentum. 
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-^orb(x/) in Eq. ((2b)) corresponds to a prograde (retro- 
grade) equatorial orbit if the final spin is aligned (anti- 
aligned). This approach agrees remarkably well with ex- 
isting simulations of moderately spinning, comparable- 
mass BBH mergers, but the assumption of mass conserva- 
tion reduces agreement with the test-particle limit to ze- 
roth order in the mass ratio q. While this approximation 
is strictly valid to 10% as noted by BKL, it artificially 
reduces xs iii the important highly spinning regime. 

To improve on BKL, we seek an analytic expression for 
the final mass Mj for arbitrary mass ratios and initial 
spins. In the spirit of their work, we propose 



A'//,k = M-m[1-S(x/)], 



(3) 



where M = ttii -I- TO2 is the total initial mass and ^ = 
Tn\milM is the reduced mass. This proposal both agrees 
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FIG. 1: The final spin parameter x/ for mergers with equal 
initial spin parameters (xi = X2) and mass ratio v. Each 
curve gives a different value of xi, which can be read off the 
y-axis as x/ = Xi for v = ^. Curves in blue (solid) denote 
our predictions, red (short-dashed) those of BKL, and green 
(long-dashed) those of HB. To-p panel: the highly spinning 
regime (xi > 0.95). Bottom panel: aligned spins (xi > 0). 



FIG. 2; The final spin parameter x/ as a function of 1/ for 
the merger of BBHs with equal and nearly maximal initial 
spin parameters (xi = X2 > 0.99). The red (short-dashed) 
curve shows the maximum initial spin (xi — 0.9916) for which 
this function monotonically increases with v. The blue (long- 
dashed) curve shows the maximum initial spin (xi — 0.9987) 
for which black holes are spun up by equal-mass mergers. 



with Eq. (jlap to first order in q and is symmetric under 
exchange of the black hole labels "1" and "2" . Dividing 
Eq. (j2bp by the square of ([3]), we find 



Xf = 
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(4) 
for the final spin parameter, where v = fi/AI. 

In Fig. [U we compare this new prediction for xf (here- 
after K) with those of HB (after including S2 in Eq. (jlbp ) 
and BKL. We consider BBHs with equal initial spins 
(Xi = X2 = Xi) aligned with the orbital angular mo- 
mentum. All three approaches reproduce the trivial ze- 
roth order result that x/ = Xi in the test-particle limit 
{ly = 0). However, only K and HB provide the correct 
first-order behavior (the value of dxf/dv) in this limit 
[2^ . The BKL assumption Mf = M implies that E =1 
and makes the denominator of Eq. ([4]) equal unity. Over- 
estimating Mf to first order in ly leads to a first-order 
underestimate in Xf = Sf/M'i, accounting for the neg- 
ative values of dxf/dv given by BKL for Xi ^ 0.95. K 
and HB predict dxj /dv{xi ^ 1) > 0) with equahty for 
maximal initial spins. This agrees with the famous result 
of Bardeen [l3| that black holes can be spun up to the 
Kerr limit x = 1 after accreting a finite mass of test par- 
ticles. Despite its success for i^ == 0, HB predicts naked 
singularities (x/ > 1) for comparable-mass mergers. 

Let's examine the predictions of K more closely in the 



limit of nearly maximal spins (xi = X2 > 0.99). In Fig.[2l 
we show X/(^) at intervals of 0.002 in Xi- For Xi ^ 0.9916 
(the red, short-dashed curve), Xfiy) is monotonically in- 
creasing implying that equal-mass mergers {v ~ 0.25) are 
the most efficient way of spinning up a black hole. As Xi 
increases, the maximum t'max of X/(^) moves to lower 
V. For Xi ^ 0.9987 (the blue, long-dashed curve) equal- 
mass mergers spin down black holes rather than spin- 
ning them up. Finally, at Xi = 1, the maximum reaches 
V ~ Q making x/(^) monotonically decreasing. Thus 
maximally spinning black holes will be spun down for all 
mergers with v > Q. The merger of maximally spinning 
equal-mass BBHs yields a black hole with x/ — 0.9988. 

This behavior is summarized in Fig.[3l At Xi ^ 0.9916, 
I'max falls below 0.25 and the red curve peals away from 
the blue curve in the top panel, indicating that equal- 
mass mergers are no longer the most efficient way of spin- 
ning up highly spinning black holes. While the blue curve 
belonging to equal-mass mergers falls below the x/ = Xi 
reference curve for Xi ^ 0.9987, the red curve always has 
Xf ^ Xi implying that BBH mergers can always spin 
up black holes provided v is low enough. Though the 
specific values provided here are only as accurate as the 
approximation itself, we believe this qualitative behavior 
will ultimately be verified by numerical simulations once 
near maximally spinning initial data can be constructed. 

Let's compare these predictions to recent numerical- 





FIG. 3: Top panel: The final spin parameter Xf ^s a func- 
tion of Xi for the merger of BBHs with equal and nearly 
maximal initial spin parameters (xi = Xa = Xi)- The blue 
(long-dashed) curve corresponds to equal-mass BBHs, the red 
(solid) curve to mergers with mass ratio I'max that maximizes 
Xf{^)j ^nd the black (short-dashed) curve to the reference 
Xf = Xi- Bottom panel: The value i-'max(Xi) described above. 



relativity simulations. Wc begin with initially non- 
spinning BBHs of mass ratio i/. The top panel of 
Fig. [5] shows the final masses Mf determined from the 
energy radiated (squares), apparent horizon (triangles), 
and quasi-normal modes (QNMs, circles) [1J|. Eq. ([3]) 
underestimates Alf by nearly a factor of two for equal- 
mass mergers, possibly indicating that as much as 50% 
of the energy is radiated after ISCO for these mergers. 
Better agreement is seen at lower mass ratios, in line with 
arguments that the energy radiated durin g m erger and 
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ringdown scales as v^ rather than v 
these poor estimates of Mf , our predicted final spins Xf 
are amazingly accurate; the equal-mass prediction {xf ~ 
0.687) lies within the numerical errors (Ax/ ~ 0.002) 
of most simulations, an order-of-magnitude improvement 
over BKL. This is telling us something profound; while 
significant energy and angular momentum are radiated 
after plunge, they are radiated in the appropriate ratio 
to preserve the spin Xf — Sf /M'i. 

This is not always the case for spinning BBHs. Fig. [5] 
shows the final masses and spins for equal-mass, equal- 
spin BBH mergers. If we believe our approximation, as 
much as 45% of the energy is radiated after the inspiral 
stage for aligned, highly spinning configurations. This 
energy loss tends to drive Xf up, but apparently angular 
momentum is radiated even more efficiently after plunge 
because the numerically determined spins are below our 



FIG. 4: Top panel: The final mass Mfiy) for initally non- 
spinning mergers. The red (short-dashed) curve is the BKL 
approximation, while the blue (solid) curve is that of Eq. ([3|. 
The square, triangle, and circular points are three numerical 
estimates of Mf [ij]. Bottom panel: The final spin Xfiy) for 
non-spinning mergers. The red (short-dashed) curve is the 
BKL prediction; the blue (solid) curve is that of Eq. @ . The 
black points are three estimates of the final spins from [l4| . 
while the green Xs are simulated final spins from [15| . 



predictions. The BKL approximation actually does bet- 
ter for these configurations, but this appears to be a con- 
spiracy between the excess final mass compensating for 
the failure to account for the angular momentum radi- 
ated after plunge. Angular momentum is radiated less 
efficiently after plunge for anti-aligned spins, leading us 
to underestimate Xf for Xi < 0- Future work will at- 
tempt to improve our approximation by accounting for 
radiation after the ISCO. 

This same trend is seen in equal-mass, unequal-spin 
simulations [18[. There is excellent agreement between 
our approach and all simulations with vanishing total 
initial spin like those shown in the upper left panel of 
Fig. [5] Our approach improves on BKL except for the 
highly spinning aligned configurations where their artifi- 
cially high Mf drives Xf down. 

BKL used conservation of angular momentum to pre- 
dict with great success the final spin in BBH mergers. 
We have improved on their approach by allowing for en- 
ergy loss during the inspiral, and have uncovered qualita- 
tively new dependence of the final spin on the spins and 
mass ratio of the initial BBHs. Equal-mass BBH mergers 
might lead to higher final spins than expected, rivaling 
the highest spins (x — 0.998) allowed by physical accre- 
tion flows |2(|. Black holes can always be spun up by 
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FIG. 5: Top panel: The final mass Mf{xi) for equal-mass, 
equal-spin mergers. The red (short-dashed) curve is the BKL 
prediction, while the blue (solid) curve is that of Eq. ([Sj. 
Square points are simulations from [101; the purple X is from 
[131 . Bottom panel: The final spin parameter x/(^) for these 
mergers. The red (short-dashed) curve is the BKL prediction, 
the blue (solid) curve is that of E g. ffl . Black squares, purple 
Xs, and green triangles are from [ly], [l7| . and |18l |. 



aligned mergers provided the mass ratio is small enough. 
We have found that for non-spinning BBH mergers, exact 
analytic results for test particles on Kerr geodesies can 
be extrapolated all the way to equal masses with better 
than 1% accuracy. Our approach provides quantitative 
predictions for regimes not yet accessible to numerical 
relativity, and helps guide the choice of simulations to ex- 
plore the 7-dimensional parameter space of BBH mergers 
formed by i'. Si, and S2 ■ This complements fitting formu- 
lae that successfully describe final spins for comparable- 
mass mergers ISl ll9|, l21| . It can help identify where in 
parameter space black holes will be spun up by mergers, 
and where the final spin will vanish leading to interest- 
ing orbital dynamics. Understanding BBH mergers is 
an essential goal of astrophysics, general relativity, and 
gravitational-wave detection, and can only be achieved 
by combining numerical relativity, post-Newtonian tech- 
niques, and analytic approximations like ours. 
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FIG. 6: Final spin Xf ^s a function of xi for the four se- 
quences of equal-mass, unequal-spin simulations presented in 
[Iq l . The value of X2 for each sequence is listed on each panel, 
and we have exchanged xi and X2 in the upper right panel 
for convenient presentation. Blue (solid) curves are our pre- 
dictions, red (dashed) curves those of BKL. 
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Correct dependence on v helps calibrate fits for x/; for 
example, Eq. @ implies the coefficients of [ly| should be 
S4 = -(2V3)/9, to = -(16y3)/9, and ii = 2^/2,. 
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